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Abstract 

It is well known that a state of pure shear has distinct sets of basis vectors or 
coordinate systems: the principal axes, in which the stress cr is diagonal, and pure 
shear bases, in which diag<x = 0. The latter is commonly taken as the definition 
of pure shear, although a state of pure shear is more generally defined by tr cr = 0. 
New results are presented that characterize all possible pure shear bases. A pair of 
vector functions are derived which generate a set of pure shear basis vectors from any 
one member of the triad. The vector functions follow from compatibility condition 
for the pure shear basis vectors, and are independent of the principal stress values. 
The complementary types of vector basis have implications for the strain energy of 
linearly elastic solids with cubic material symmetry: for a given state of stress or 
strain, the strain energy achieves its extreme values when the material cube axes 
are aligned with principal axes of stress or with a pure shear basis. This implies 
that the optimal orientation for a given state of stress is with one or the other vector 
basis, depending as the stress is to be minimized or maximized, which involves the 
sign of one material parameter. 

1 Introduction 

A state of pure shear is characterized by 

tr<r = 0, (1) 

where cr is the (symmetric) stress tensor. The principal axis form of stress is 

a = 0i ei <g) ei + a 2 e 2 <g> e 2 + a 3 e 3 <g> e 3 , (2) 

where the orthonormal triad {e 1; e 2 , e 3 } is unique up to reflections, i.e. {±e 1; e 2 , e 3 }, etc. 
In a state of pure shear the three principal stresses satisfy 

ci + o"2 + 0-3 = 0. (3) 

The condition ((TJ) is unaltered under a change of basis, say from {ei, e 2 , e 3 } to {a, b, c}. 
We define a pure shear basis as one for which the diagonal elements of stress vanish, 

a aa = a ■ (era) = 0, cr bb = b • (crb) = 0, a cc = c • (crc) = 0. (4) 
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Coordinate systems with diagcr = have been the subject of recent interest [U 121 El, 
focusing on the relation between eqs. J!} and (J3J), and on the characterization of all 
possible pure shear bases. The three conditions Q are the conventional definition of pure 
shear, although Q is a more general starting point. The equivalence is well known, and 
is mentioned, for instance in Love's book (Section 16 of [I]). Perhaps the first proof that 
eq. implies eq. (j3J) was by Gurtin j3] , who in turn ascribed the method to Lew. Belik 
and Fosdick P offered two distinct proofs, geometric and algebraic. They showed that 
there is a one parameter set of pure shear coordinate systems. Boulanger and Hayes [2] 
subsequently gave an interpretation in terms of elliptical sections of the stress ellipsoid. 
Ting |B] provided a characterization of directions n such that <j nn = in terms of the 
total shear in the plane normal to n, r = |trn|. 

The first objective here is to provide a simple means to characterize all possible pure 
shear bases. What distinguishes the present discussion from the recent notes 121 E] is 
that the pure shear basis vectors are determined in a manner that is, as far as possible, 
independent of the values of the principal stresses, oi, o 2 and o 3 . We show that one 
member of the orthonormal triad {a, b, c} determines the other two through a cyclic 
vector function defined relative to the principal stress axes. This representation is achieved 
by considering the conditions necessary for a given orthonormal triad to be a pure shear 
triad. The conditions are independent of the principal stresses, and the vector triad is 
completely defined in terms of one element, a, say. The stress then enters through the 
condition that o aa vanish, which is satisfied if a lies on the intersection of an elliptical 
cone and the unit sphere. In this formulation no assumptions, such as the sign or ordering 
of the principal stresses, are required. The complete set of pure shear bases follows simply 
from the vector functions which define the triad. 

The principal axis and pure shear coordinate systems for a state of pure shear have 
direct implications for strain energy in materials with cubic symmetry. Thus, the extreme 
values of the strain energy correspond to aligning the cube axes with the two distinct sets 
of orthonormal base vectors: the principal stress axes, and the pure shear axes. This 
provides a simpler method than those discussed in jH] for characterizing the orientations 
of the cube axes which maximize or minimize strain energy for a given state of stress. 

We begin in Section [2] with the general result for the pure shear orthonormal bases, 
which is discussed subsequently. Implications for strain energy are explored in Sectional 

2 Pure shear axes 

2.1 Generic functions for the basis vectors 

Consider basis vectors {a, b, c} such that (@J) is satisfied, i.e. 

a\a\ + a 2 a\ + a^al = 0, (5a) 
<j x b\ + a 2 b\ + 0-363 = 0, (5b) 
a\c\ + a 2 c\ + (X3C3 = 0, (5c) 

where components relative to the orthonormal principal axes 

of stress satisfying [ei, e 2 , e 3 ] = 1, i.e. a = a±ei + a 2 e 2 + a 3 e 3 . Equations (J5J) imply that 
each of the unit vectors lies on the curve C defined by the intersection of an elliptical cone 
with the unit sphere: 

C = {xiei + x 2 e 2 + x 3 e 3 } : a x x\ + a 2 x\ + 0-3X3 = 0' + £3 + £3 = 1. (6) 
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Alternatively, writing the three conditions as a matrix equation, 



b 2 

"2 




(7) 



it is clear that the following is a necessary condition for the triad {a, b, c} to simultane- 
ously lie on C, 



4 



h(&, b, c) = det b\ b 



0. 



(8) 



We will show that if one member of the triad is known then this compatibility condition 
defines the other two. 

The main result is as follows: 

Theorem 1 The triad {n, v^n), v 2 (n)} is orthonormal with [n, v 1; v 2 ] = 1 and satisfies 
the compatibility condition for shear axes, i.e. h(n, v 1; v 2 ) = (eq. $B)). The vector 
functions v x and v 2 are defined 



n 



^(u_-su + ), v 2 (n) = ^(u 



v^2 



S U i 



(9) 



where 



s(n) = sgn [nin 2 n 3 (nl - §)(n 2 - ~){n\ - §)] , 
and the orthogonal unit vectors u ± are 

p±ni p±n 2 p±n 3 



u±in) 



P±l n J 



n 2 3 ±g 



ei + 



+ 



n\ ±g 



e 2 + 



+ 



n\ ±g 



e 3 , 



117, 



n- 



±g) (nl-n\±g) (n\-nl±g)' 



-1/2 



g(n) = [l- A(n\nl + n\n\ + n\nf) + §n\n\n 2 z \ 



(10a) 

(10b) 
(10c) 



The vector functions Vx and v 2 together with n form an orthonormal triad, and there- 
fore satisfy the cyclic identities 



vi(n) = v 2 (n) x n, 



v 2 (n) = n x vWn) 



n = vi(n) x v 2 (n). 



An immediate corollary of Theorem 1 is that the triad {n, m, p} form a pure shear 
basis if one lies on C and the other two are given by Theorem^ For instance, if n satisfies 



ain\ + a 2 n\ + a 3 n\ 



0. 



and m and p are defined by the vector functions, 

m = vi(n), p = v 2 (n) 



(12) 
(13) 



then o~ nn — o~ mm — o~pp — 0. 

Proof of Theorem n Consider eq. (jSJ) for an orthonormal triad {n, m, p}. Since the 
vectors are orthonormal it follows that the sum of the three elements in each column and 
each row of the determinant is unity. Thus, 



/(m) 











n\ 


n\ 


2 






m\ 




ml 


1 


1 


1 



I 2 2\ 2 i / 2 2\ 2 i / 2 2\ 

[n 3 - n 2 )m 1 + (n x - n 3 )m 2 + (n 2 - njm 



(14) 
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is a necessary but not sufficient condition for the transformation from principal axes to 
pure shear axes. We will now use this condition to find m and p for a given n. Consider 
the generalized function 

F(u) = /(u) - Auu + 2pnu. (15) 

The form of F is motivated by the compatibility condition for pure shear axes augmented 
by terms that constrain the vector u to be of unit magnitude and orthogonal to n. The 
Lagrange multipliers A and p are determined by requiring that / is stationary with respect 
to u for a given n. Setting the gradient to zero implies 

pni pn 2 pn 3 

u= — TTT ei + ^ TTT e i + — TTT e3 ' ( 16 ) 

n 2 ~ n i + ^ nl — n{ + X n{ — n 2 + A 

where p is such that u is of unit magnitude, and A follows from the condition n ■ u = 0, 

U " ' U " n * -0. (17) 



n\ — n1 + \ n| — n\ + A n\ — n 2 + A 
This is a quadratic equation in A which can be simplified using |n| = 1 to the form 

A 2 = 1 - ^{n\n\ + n\n\ + n\n\) + §n\n\n\. (1 



Thus A = ±g(n), where g is defined in eq. (jlOcj) . 

Equations (|15|) through (fTTj) imply that the function F(u) vanishes at the stationary 
value of u, and hence 

/(u) = A. (19) 

The function /(u) therefore achieves its maximum and minimum values ±g(n) along the 
directions vectors u± defined by (HOaj) . It may also be checked, again using (fTTj) . that 
u + and u_ together with n form an orthonormal triad. Furthermore, since f(m) is a 
quadratic form in the vector m, it must have the form 

/( m ) = 0(n)[(m- u + ) 2 - (m- u_) 2 ] 

= — 2g(n) m • Vi(n) m ■ v 2 (n) for m ■ n = 0, (20) 

where Vi and v 2 are the orthogonal bisectors of the directions u + and u_, and therefore 
correspond to the directions for which / vanishes. This proves the Theorem. 

The parameter s — ±1 is introduced to ensure that the scalar triple product [n, v 1; v 2 ] = 
n • (vi x v 2 ) is +1. This guarantees that the transformation from principal axes to the 
pure shear basis is a proper orthogonal transformation. 



2.2 Discussion 

2.2.1 Stresses dual to the principal stresses 

Before discussing the vector functions in Theorem 1, we note some alternative formu- 
lations. We first introduce the stress parameter ov^) for the pure shear axis n. Thus, 
using the pure shear definition (jH)), the condition (fT2j) for the vanishing stress a nn may be 
written 
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which defines <7( n ). Equivalently, 



ni 




™1 = 5 +77— ^1-^3 . (22) 



The stress parameter <T( n ) can be related to the total traction r n , defined by r 2 = \crn\ 2 
or using a nn = 0, 

2 2 2, 2 2 , 2 2 / no \ 

r„ = o x n x + a 2 n 2 + a 3 n 3 . (23) 

Thus, 



J_ = 3r w 2 - (a 2 + o\ + gj) 

O-(n) - CT2)(0"2 - cr 3 )(cr 3 - 0"i) ' 

and conversely 

^ = o^i 2 + CT 2 + ^) + ~ °2)((T2 ~ cr 3 )(a 3 - m). (25) 

O OCT( n ) 

Note that elimination of at n ) from eq. (|22|) gives the known relations for pure shear [3] 

n? - T ' + ^ 



1 (ffl - 02){O1 ~ ^3)' 

"1 = . "" + ,, g3 " (26) 

(cr 2 - 0- 3 )(cr 2 - 0"l) 
^2 = T n + glgl 

3 (er 3 - 01) (cr 3 - cr 2 )' 

Let {a, b, c} be a pure shear basis satisfying eq. (HJ), and define <7( a ), ovm and <7( c ) 
according to eq. (}2~T]) . It follows from eq. (J2"3j) that 

1 , 1 , 1 _3 (r a 2 + T- 2 + r 2 -a 2 -a 2 -a 2 ) ^ (2?) 



C"(a) CT(6) CT(c) — 0"2)(0"2 — ^3)(C3 — 

However, t 2 + r 2 + r 2 = tr ([a ® a + b ® b + c ® c]cr 2 ) = tr cr 2 , since {a, b, c} is an 
orthonormal basis. We therefore deduce the identity: 



— — I — - — I — — = for a pure shear basis. 

CT (a) °"(b) °"(c) 



(28) 



This condition for a pure shear basis should be compared with the condition (JHJ) satisfied 
by the principal stresses. The stress parameters <T( a ), <T(j,) and cr( c ) are in this sense dual to 
the three principal stresses, and the condition (J2"%|) is the dual of the pure shear condition 
tr cr = 0. 

2.2.2 An alternative form of the vector functions 

Returning to the proof of Theorem 1, specifically the expression (fTH|) for the vector u at a 
stationary value of the functional, we note that terms such as (n 2 — nf) may be replaced 
by (T\/(T(n) using eq. (|21j). As a result, stationary values of the vector u are of the form 

Pm P n 2 P n 3 (on , 
u = — -ei + ~^ e 2 + — re 3 , 29 

O"! + A cr 2 + A cr 3 + A 
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where A is now determined by the necessary condition n • u = 0, 



+ 



0. 



(7i + A (7 2 + A a 3 + A 
This can be simplified using the fact that a nn = to 

A 2 + a 2 a 3 n\ + a^^l + Oia 2 n\ = 0. 
Consequently, the intermediate vectors u± in Theorem 1 can be expressed 



u±(nj 



d(n) 



0"! ± d 



ei 



2 i ^ 



e 2 



cr 3 ± d 



e 3 , 



II', 



+ 



rio 



Y(a 1 ±dy (a 2 ±dy (a 3 ±dy 

1/2 



-1/2 



[nl<j 2 <J 3 + ^2°"30-l 



(30) 

(31) 

(32a) 
(32b) 
(32c) 



These formulae provide an alternative representation which combines the vector coordi- 
nates of n with the principal stress values. We are now in a position to discuss the vector 
functions Vi and v 2 . 

2.2.3 The vector functions 




Figure 1 : All directions are covered by considering n restricted to the irreducible 1 / 48th of the 
cube surface shown. The region is defined by the triangle on the cube surface with vertices in 
the 001, 110 and 111 directions. 



The vector functions Vi and v 2 of Theorem 1 are expressed in two different forms 
through eqs. (jlOj) and (}3*2*|) in terms of the secondary vector functions u_|-. It is useful to 
keep the alternative representations in mind as we consider the dependence of the vector 
functions on n and on the principal stresses. The range of variation of n may be reduced, 
with no loss in generality, to the interior of the triangle depicted in Figure d The vertices 
001, 110 and 111, described using crystallographic notation, correspond to the generic 
directions e^, (e$ + ej)/y/2 and (ei + e 2 + e 3 )/\/3 (i ^ j =1, 2 or 3). 
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The functions vi(n) and V2(n) of Theorem 1 are well defined and unambiguous for n 
lying within the interior of the triangle although some care is required when evaluating 
the functions on the boundary. Consider, for instance the boundary between the 001 and 

110 vertices, which is characterised by one of the rectangular components of n tending to 
zero. Consider n 3 =o(l) with (n\ — n%) =0(1), then to leading order in n 3 the vectors u± 
are e 3 and e 3 x n. This approximation holds for all n along the edge between 001 and 110 
in Fig. ^ including n approaching a principal axis (the 001 direction). As n approaches 
a face diagonal (the 110 direction) the parameter g — > and consequently the function / 
of eq. (|2()jl is identically zero for all vectors u in the plane orthogonal to n. In this case, 
the vectors u± and hence vi and V2 are degenerate. Note that both the 001 and the 110 
directions are pure shear axes if and only if one of the principal stresses vanishes. Thus, 
if n = e 3 is a pure shear axis then a nn = 03 is zero. Similarly, if n = (ei + 62) / v^2 is 
a pure shear axis then o\ + 02 = 0, which implies 03 = since trcr = 0. We note that 
< g < 1, with g(ei) = 1, g{^e t + ^e,-) = for i ^ j, and 0(^1 + ^e 2 + ^e 3 ) = 0. 
In general, the function g is zero only along face and cube diagonals. 

The 111 direction is a unique and special direction in terms of pure shear axes. First, 
we note that regardless of the principal stresses, 111 is always a pure shear axis by 
virtue of the constraint (jH)). At the same time, the function g(n) vanishes for n in the 

111 direction, and the formula (jlO)) is not defined for this particular n. However, the 
alternative formulation of eq. (J32j) yields 



u±(lll) = ai±a ° a2±ao a3±ao , (33) 



1 . 1 + 1 



(oiio-o) 2 ^ (<X2±<7 ) 2 ^ (o-3±(T0) 2 

1/2 



where the positive stress cr = <i(lll) = [ — (<7 2 cr 3 + cr 3 o"i + <7icr 2 )/3] may be written 
a = -j=t with t q = [-(<t? + o\ + 0$)] 1/2 . (34) 

t is the total traction for the 111 direction (r = r n ). In summary, (u + (lll)±u_(lll))v / 2 
together with 111 define the pure shear basis. 
Using the general relation 

»(n) = ^, (35) 

it follows that 1/V( n ) vanishes in the 111 direction: 

1 



0. (36) 



This is also apparent from the relation (|2~2~|) at n = 111. Equation (|22jl. which can 
be considered the definition of the dual stress cri n \, also provides a parameterization of 
possible directions n in terms of <7( n ), considered as a variable. The special direction 111 
corresponds to 07 re ) — > 00, and neighbouring n directions can be found by continuation. 
By starting at the 111 direction one could develop all pure shear bases for a given state 
of stress. Thus, based on eq. let 



1/ ^ + 3P1 ei + d ^ + -P2 e 2 + d ~ + ^P3 e 3 , (37) 
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where the unit vector p is defined by the principal stresses 

^ (gg ~ gg) o , (gi-ga) ^ , (o-2-o-i) „ , QS ^ 

p = ei H e 2 H e 3 . (38) 

t"o To t 

The range of the parameter t is such that the arguments of the square roots lie between 
and 1. The complete set of pure shear bases for a given state of stress is then characterised 
by the set of triads {n(t), vi(n), v 2 (n)} defined by all permissible t. 



3 Strain energy in cubic materials 

We now consider the relationship between the distinct types of basis for a state of pure 
shear and the strain energy in a solid. The strain energy U is assumed to be quadratic 
in the symmetric strain tensor e, but using the linear relation between stress and strain 
the energy can be written as a quadratic in stress. Our focus is the dependence of the 
strain energy on the orientation of the principal axes of the stress <r, which is no longer 
constrained to be a state of pure shear. An alternative point of view is that for a given 
state of stress or strain, we consider how the strain energy varies as the crystal axes are 
rotated. The general question of how strain energy behaves in this manner and what 
are the optimal crystal orientations has been considered in several papers. Rovati and 
Taliercio [2j provide a good review of the topic, and they also examine the particular 
case of materials with cubic symmetry. The author considered the same question and 
derived several results concerning orientations of the cube axes that yield extremal values 
of strain energy. The purpose of this section is to revisit this issue afresh, and show that 
the question of which orientations provide the extreme values of U is easily and directly 
answered using the notion of pure shear bases. 

The isotropic invariants of strain are Ii SO (e) = tre and II\ so (e) = tre 2 , and invariants 
with cubic symmetry are / C ub(^) = tr = Ii SO (e) and J/ cu b(s) = tr(De) 2 . Here, D is a 
fourth order tensor defined by the cube axes, which are assumed to be {a, b, c}, 

= a<g>a<g>a<g>a+b<g>b<g>b<g>b + c<g>c<g>c<g>c. (39) 

The strain energy in the linear elastic theory is 

2U = Kll(e) + 2/i! [J/ iso (e) - // cub (e)] + 2^ 2 [ll cuh (e) - ^/ 2 so (s)] . (40) 

This is positive definite if and only if the three constants k, \L\ and /i 2 are all positive, 
which is assumed. These parameters, called the "principal elasticities" by Kelvin [H], 
are related to the standard (Voigt) notation for stiffness: K = (en + 2c i2 )/3, /ii = c 44 , 
A*2 = (en — Ci 2 )/2. The fourth order tensors of elastic stiffness and its inverse, compliance, 
follow from above. Thus, 

a = Ce, £ = So-, (41) 

where S = C -1 and 

C ±x = (3k) ±1 J + (2/i 1 ) ±1 L + (2/i 2 ) ±1 M. (42) 
The tensors are as follows: JS = |(trS)I for all S G Sym, 

L = I - D, M = B - J, (43) 
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and I is the fourth order identity, IS = S for all S G Sym, which can be written 

I = J + L + M. (44) 

As a consequence, the strain energy can be written in an alternative form as a quadratic 
function of the stress invariants: 

2U = ^U°) + ^ [^isoW - I/ Cub (<r)] + ^ [/JcubW - \lL{<r)] ■ («) 

The stress may be partitioned into hydrostatic and deviatoric parts, 

- = ^/ iso (o-)I + o- / , (46) 

where the deviatoric stress satisfies tr cr' = 0, and hence can be considered as a pure shear. 
Using /cub(f) = -Aso(c)5 it follows that 

U iso (<r) - ij2o(<r) = // iso (<r'), ^cub(^) - ii2o(«r) = U C ub(^). (47) 



Hence 



217 = ^°< CT > + w, Ih ^ )+{ w, - w) 11 ^' (48) 

and all the orientation dependence is through the final term. Note that 

< H cuh (*') = a'J + a' bb 2 + a'J < II iB0 (<r') = a[ 2 + a' 2 2 + a' 2 , (49) 

where a' aa = a • (tr'a) etc., and a[ are the principal deviatoric stresses, satisfying a[ + 
°2 + °3 = 0- The maximum value of II cn ^(cr') is obtained when the cube axes {a, b, c} 
coincides with the principal axes {ei,e 2 ,e 3 }. The minimum value i7 cu b(<r') = occurs 
when the cube axes {a, b, c} form a pure shear basis. The extreme values of the strain 
energy are therefore 

f 2/7T' {a, b, c} = pure shear axes, 
2t/ = -/?»+JJ iso (tr0x I W (50) 



9k 



j-j, {a, b, c} = principal axes. 



These results are unchanged by transformations congruent with cubic symmetry, i.e. 
{a, b, c} — > {—a, b, c}, etc. Hence we have the second main result: 

Theorem 2 If \i<x > \i\ (fi 2 < then the strain energy is maximum (minimum) when 
the cube axes {a, b, c} are aligned with a pure shear basis for the deviatoric stress. The 
strain energy is minimum (maximum) when the cube axes coincide with the principal axes 
of stress. 

These results may be phrased in terms of strain rather than stress. 

The state of pure shear may be reconsidered in terms of the cubic fourth order tensors. 
Thus, a state of pure shear, tr cr = 0, is equivalent to 

Jcr = 0, Ker = cr, <^ pure shear, (51) 
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where the fourth order tensor K is 

K=I-J=L+M. (52) 

We note in passing that a state of pure hydrostatic stress (pressure) is the dual of eq. 
(l5Tj) : Jcr = cr, Kcr = 0. Concentrating on pure shear states of stress, the second partition 
of eq. (|52j) enables us to define the distinguished basis vectors as follows 

Lcr = cr, Mcr = 0, <^ {a, b, c} = pure shear axes, (53a) 
Lcr = 0, Mcr = cr, <^ {a, b, c} = principal axes. (53b) 

The fourth order tensors J, K, L and M act on elements of Sym, which is a 6-dimensional 
space. In this way, dim J = 1, dim L = 3, dim M = 2 and dim IK = dim L + dim L = 
5. Thus, a state of pure shear, (|HT|) . is characterized by 5 degrees of freedom, while 
hydrostatic stress has only a single degree of freedom. For the special and complementary 
pure shear configurations of (|5*3*|) . the principal axes state is defined by two parameters, 
since <Ti, cr 2 and <t 3 are not independent but satisfy the pure shear constraint eq. (JSJ). 
The pure shear axes have an extra degree of freedom, which is reflected in the non-unique 
nature of pure shear basis, and in the fact that they define a one parameter family of 
coordinate systems, as discussed by Belik and Fosdick £Q. 

Finally, the relationship between the deviatoric stress and strain is simple at the energy 
extrema. Thus, eqs. ()42)1 and (}53|) imply that cr' = 2fiie' if the cube axes are pure shear 
axes, and cr' = 2^2^' when they are the principal axes of stress and/or strain. 
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